Phase Diagram of the S =~ Frustrated Coupled Ladder System 
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We present a theoretical study of the magnetic phase diagram of the frustrated coupled ladder 
structure realized recently in several materials. This system displays a nondegenerate spin-gap state 
in the dimer limit and an infinitely degenerate spin-gap state in the regime of weakly-coupled zig- 
zag chains. Between these we demonstrate the existence of gapless, magnetically ordered regions 
whose order is antiferromagnetic close to the honeycomb lattice limit, and incommensurate along 
the chains when all three magnetic interactions compete. 

PACS numbers: 75.10.Jm, 75.30. Kz, 75.40.Cx 
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Spurred by developments in the field of high- 
temperature superconductivity, rapid progress is now be- 
ing made in the preparation of materials with similar at- 
tributes, and the potential for fascinating new physics. 
These are inorganic, low-dimensional quantum magnets, 
where the relative strengths of the magnetic interac- 
tions, which in most cases give antiferromagnetic cou- 
pling between ions with spin S = i, result in systems 
such as chains and ladders, which are effectively one- 
dimensional (Id), or 2d coupled ladders and depleted 
planes. One structure of particular interest is the frus- 
trated coupled ladder, shown in Fig. 1. This conforma- 
tion is realized in the cases of the prototypical "ladder" 
compoundEI SrCu203, where the frustrating interaction 
Ji is ferromagnetic (FM) and small, so that the J2-J2 
ladders are only weakly coupled, and in the "zig-zag" (or 
J1-J2) chain compoundEI SrCu02, where the interchain, 
coupling J 2 is small. The depleted planar compoundEI 
CaV205 represents a system where all three interactions 
are thought to be of similar strength. 
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FIG. 1. Schematic representation of the frustrated coupled 
ladder system. 

The phase diagram obtained on alteration of the pa- 
rameters Ji , Ji and J 2 is particularly rich. The two-chain 
ladder obtained at J\ = has a singlet ground state of 
resonating dimers with a gap to all spin excitations, while 
the zig-zag chain at J 2 = shows for J 2 > Jic a doubly 
degenerate ground state of alternating dimers which also 
has a gap. What appears between these limits remains 
poorly understood, and is one subject of the current anal- 



ysis, but we may gain an initial indication from the hon- 
eycomb lattice obtained when J 2 = 0, where the gapped 
J 2 dimer phase is replaced by a magnetically ordered 
state at a critical value of increasing J\. We investigate 
the system using a variety of analytical and numerical 
methods, proceeding in (a) from the limit of ladder rung 
dimers, in (b) from that of weakly coupled zig-zag chains 
and in (c) by considering the ordered phases which are 
found to occur between these limits. 

a. Dimer Limit: in the limit of large J 2 , the system 
has a non-degenerate ground state, whose wavefunction 
is a product of rung dimer singlets, with an energy gap 
to spin excitations. This state is well spited to exami- 
nation by the bond-operator technique ,u which is based 
on transforming the four spin states on each rung to one 
singlet and three triplets. The gap between these gives 
the stability of dimer order. The method has also been 
showna to be applicable to spin ladders (vanishing J\) 
for values of J2 up to J 2 , and can be expected to have 
similar validity in the case of the anisotropic honeycomb 
lattice (vanishing J 2 ). 

Following the treatment of Ref. || for a unit cell con- 
taining two dimers, and retaining in the spin-gap regime 
terms only to quadratic order in the triplet spin exci- 
tations t\ a , we obtain a system of threefold degenerate 
magnons with energy spectrum 



1 + dai 



where d = 2J 2 s 2 / (\.J' 2 - /i), 



= A cos k z ± A' cos ik x cos hk z , 



(1) 



(2) 



s denotes the magnitude of the singlet condensate, /i is 
the global chemical potential, and we define A = J2/J2 
and A' = J\j J 2 . Solution of the mean- field equations at 
T = gives the most important property characterizing 
the non-degenerate ground state, the spin gap 



A = (i J' 2 - p) y/l-df(\,\>). 



(3) 
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For Ji > 4J 2 , / = A' — A = a k min at the commen- 
surate wave vector kyi = (0,0), while for J 2 > \J\ 
one has / = A + ^ occurring at an incommensurate 

= (0,2cos _1 jr). In the spin-gap regime, the maxi- 
mum of the static structure factor S(k) coincides with the 
minimum of the gap, and can be shown to move to incom- 
mensurate wave vectors with increasing J 2 beyond \J\- 
The boundaries of the spin-gap phase are found where the 
A — > 0, and are shown in Fig. 2. When J\ > 4J 2 , from 
kyi one may deduce a FM arrangement of rung spin pairs 
(| TD), which corresponds to antiferromagnetic (AF) or- 
der along the J\ chains. For J 2 > jJi, describes spi- 
ral order along the chains with exactly the pitch known 
from the classical solution for the J\- J 2 chain. The phase 
boundaries in Fig. 2 appear well outside the regime of 
validity of the bond-operator technique, which overesti- 
mates the stability of the dimcr state. However, they 
remain qualitatively correct, particularly in showing the 
instability of the dimer liquid to the different ordered 
phases, and furthermore provide a useful indication of 
where to apply alternative approaches. 
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FIG. 2. Phase diagram of the coupled dimer system, shown 
in the plane of J2/J1 and J2/J1, as deduced from the 
bond-operator method. 

b. Frustrated Chain Limit: for large J\ or J 2 , the 
system consists of weakly coupled chains, whose prop- 
erties are well known from the Bethe Ansatz solution. 
However, how these are affected by the weak couplings, 
which may or may not be frustrating, is an extremely 
subtle issue to which we shall return below. Here we 
consider the limit of small J' 2 , where the problem is that 
of the J1-J2 chain. This is a classic example of frustra- 
tion in Id, and we review briefly its properties (see Ref. 
|, and Refs. (12)-(18) therein). At the Majumdar-Ghosh 
(MG) point Jd = | Ji, there is an exact dimer wave func- 
tion solution!] with a twofold degenerate ground state 
and gapped excitation spectrum.cl The wave function is 
composed of frozen singlets on alternating J\ bonds, and 
has two possible realizations differing by a translation of 



one unit along the.-|Chain. The gap decreases exponen- 
tially with J 2 — J 2c EI on approach to the conformal point 
Jic/Ji — 0.2412, below which the frustrated chain re- 
tains the gapless spectrum of the nearest-neighbor spin 
chain. For J 2 > \,J\ the gap increases to a maximum 
at 0.6 < J2/J1 < 0.7, then tends towards an exponen- 
tial decay with — J 2 / Ji foxdarge J 2 . The static spin-spin 
correlation function S(k)t3 has a maximum at k = it for 
Ji < \J\i an d for J 2 > \J\ is maximal between k = ^7r 
and 7r, indicating predominant spiral correlations. 

We examine the lifting of degeneracy and the clos- 
ing of the singlet-triplet gap A by an adapted pertur- 
bational approach, which due to the degenerate nature 
of the ground state is in a strict sense variational. The 
Hamiltonian is written as H = Hq + H' , where 



Hq — J\ ^ + Ji ^ SjjSj +2i f (4) 

H' = J 2 §i,iSj,i+i + J 2 2_j ^U^iV+i! (5) 



i,l even 



i,l odd 



and in H' both i and / take either even or odd values 
only. We construct the variational wave function 



= \S)o (8 ... (8 |S%-i ® \ t )j ® 



+1 



(6) 



in which \S)i — r)\So) + r/\S„) is a linear combination 
of the two degenerate ground states on the z-th chain. 
\ T )j = J2 a (£.k,a\T k , a } + £fc+7r,a|7fc+7r,a)), where \T k>a ) is 
the eigenstate with momentum k of the isolated zig-zag 
chain, denotes a triplet excitation of the j-th chain. In 
the thermodynamic limit these form a continuum above 
the gap, and their derealization by H' gives a kinetic 
energy gain which reduces A. 

In the regime Ji <\J\ where the spin correlations are 
predominantly antiferromagnetic, the lowest-lying triplet 
excitations appear at k = and k = it. Here 77 or rf = 0, 
and the "gap equation" for J 2 takes the simplified form 
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Ji 
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where es and e Q are the ground and excited state en- 
ergies, and may be evaluated in terms of the spin-spin 
correlation function 

ImS? (w.fc) =Y,\(Tk +q ,a\S* k \S q )\ 2 5(uj-u k+q ), (8) 

a 

in which tOk+q = £k+q.a — f'S- These expressions are cal- 
culated using a Lanczos algorithm for chain lengths up 
to N — 22 sites. Setting A — yields the value J' 2c 
where the ordered phases appear. J 2c scales with the un- 
perturbed gap magnitude Ao, and vanishes at J 2 = J 2c - 
For J 2 > hJi, the situation is more complicated because 
the lowest triplet excitations appear at intermediate mo- 
menta, corresponding to predominant incommensurate 
spin correlations.Ej There is no decoupling as above, ne- 
cessitating solution of the full J 2 equation and calcula- 
tion of the energy by minimizing with respect to r\ and 
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rj . A — > for a value of J' 2 whose dependence on Ao 
varies because of the incommensuration, as shown in the 
inset of Fig. 3. That gap closure occurs at a wave vector 
away from n we take as an indication of spiral magnetic 
order. The method of solution is most accurate where 
the gap of the dimerized phase is large, and we display 
in Fig. 3 the phase diagram for the parameter regime 
0.4 < J%fJ\ < 1. The nature of the magnetic order sta- 
bilized by increasing _Zp corresponds to the wave vector 
kui maximizing S(k)t3 



0.6 



0.4 



0.2 



Spiral LRO 



Spin gap 
Dimerized Chains 



Neel LRO 



-> 4 



0.3 0.5 0.7 0.9 



0.5 



1.5 



Jo/J 



2' d 1 



FIG. 3. Phase diagram of the coupled zig-zag chain sys- 
tem, as deduced from the variational method. Solid symbols 
denote the commensurate solution (J2 < §Ji) and empty 
symbols the incommensurate. 

We may analyze the energy shift and splitting of the 
ground state by applying second-order degenerate per- 
turbation theory. At this level one need only consider 
two coupled chains, which have a four-fold degenerate 
ground state, and the energy correction per site is a 4 x 4 
matrix conveniently calculated using the Lanczos algo- 
rithm. The correction is found to be proportional to the 
unit matrix, so the degeneracy is not lifted in the thermo- 
dynamic limit to second order. This result suggests that 
the degeneracy of the ground state will be reflected in the 
presence of low-lying singlet states in the singlet-triplet 
gap for small J' 2 . We note that this behavior contrasts 
with the lifting of degeneracy which occurs in a system 
of gapped spin chains where i and I may take all values 
in Eq. (§. 

c. Ordered Phases: these may be treated by ex- 
panding in fluctuations around a state with fixed mo- 
ments, of periodically varying orientation, on each site. 
We present here results from the Schwinger boson and 
linear spin waya jricthods. In the Schwinger boson 
transformation,EM3 ordered solutions are described by 
a Bose condensation where (b) ^ becomes a parameter 
in the system of mean-field equations, and the chemical 
potential [i is determined by the requirement that the 
excitation spectrum is gapless (Goldstone modes). We 
find both long-range order of Neel type, where all cor- 
relations are AF only, and spiral order where there are 
one condensate and five bond-order parameters in the 



solution. The latter number arises because all links J\ 
are equivalent, as are all links J 2 , i.e. no solutions can 
be found where there is a symmetry-breaking analogous 
to the dimerization in the MG state, so in general only 
2 order parameters are required for each type; J 2 links 
are AF only, requiring one. Disordered solutions have no 
Bose condensation, a gapped excitation spectrum and a 
six-parameter solution including fi. The phase diagram 
is shown in Fig. 4. It is gratifying to find that the AF 
and spiral ordered phases appear in the expected loca- 
tions, and that the phase boundaries are in qualitative 
agreement with those deduced from both the dimer and 
zig-zag chain limits. We note in particular the two fea- 
tures that i) there is an upper boundary for the regime 
of spiral order at sufficiently large Ji and ii) there is a 
continuous transition directly from AF- to spiral-ordered 
regimes. 

Applying the linear spin wave technique, again (Fig. 4) 
we find regimes with both AF and spiral long-range or- 
der, and again there is a definite upper limit in J 2 on the 
latter. In this case there is no direct transition between 
the ordered phases, which are always separated by a dis- 
ordered region. This can be shown in detail by studying 
the ground state energy and staggered magnetization. 
A similar contradiction has been found previouslylia in 
the square lattice with frustrating next-neighbor superex- 
change. For this feature we have no way of distinguish- 
ing between the phase diagrams of Fig. 4 on the basis 
of these studies alone, and leave as open the possibility 
that the transition between ordered phases may be con- 
tinuous along a finite line in parameter space. We may 
not exclude the existence of a phase which is none of 
the four discussed hitherto, as an intermediate between 
both the two ordered phases and the two (degenerate and 
non-degenerate) gapped phases. 
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FIG. 4. Phase boundaries of the ordered regions, deduced 
by the Schwinger boson and linear spin wave techniques. 

All of the approximations applied have tended to over- 
estimate the extent of the parent region, although the 
deduced phase boundaries are in good qualitative agree- 
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ment. The quantitative aspects may be addressed by 
numerical techniques. When J 2 is zero, the system is 
an anisotropic honeycomb lattice, which has an unfrus- 
trated ordered state for suitably large J\ (J 2 / 0), and 
a spin-gap phase otherwise. Because the absence of frus- 
tration removes the sign problem, the quantum critical 
point on this axis can be found essentially exactly by 
large-system Quantum Monte Carlo studies. The re- 
sult J 2 / ' J\ = 1.74tL3 fixes the unknown end of the phase 
boundary of the non-degenerate spin-gap regime. Com- 
parison with the bond-operator (1.15) and Schwinger bo- 
son (2.65) results indicates the validity of each. We have 
in addition performed exact diagonalization (ED) stud- 
ies on small systems, from which for the present purposes 
it is possible to locate crossover regions between phases 
by analyzing the quantum numbers of the ground states. 
This work remains in progress, and the results will be 
presented in detail elsewhere. 

We conclude with a more schematic discussion of the 
important limits of large J\ and large Ji- When J\ is 
large, the system is one of weakly coupled chains with a 
small, marginally relevant next-neighbor frustration J 2 . 
The action of the interchain coupling J 2 has been dis- 
cussed extensively, m, recent literature for the anisotropic 
square lattice.li3li.llLa Despite the difference in behavior 
noted above for this geometry when J2 > J2c, we may, 
following Ref. [l4|, cast the interchain interaction as an 
effective, staggered mean field and deduce that any finite 
J 2 leads to magnetic order while the spectrum of the iso- 
lated chains remains gapless. Thus for any frustration 
Ji < J2c, the presence of a J 2 term can be expected by 
this argument to induce AF order. It remains to com- 
pute the value of J-ic on a frustrated chain in the presence 
of the effective staggered field: if this quantity does not 
track the phase boundary of the multiply degenerate MG- 
type ground state, we are presented with the existence 
of an intermediate phase in the region of the conformal 
point and the known phase boundary. This may be the 
state of no magnetic-order and no spin gap arising in the 
above discussions,ll§E3 which, following the analysis of 
the ordered states, may also be a candidate intermediate 
regime between AF and spiral order. 

At large J2 the problem is one of chains with compet- 
ing ladder and frustrated couplings. The ladder coupling 
is relevant, and opens a gap A linear in J' 2 j ' J2 to the 
nondegenerate state. The zig-zag chain coupling opens 
an exponentially small gap3 A ~ exp(— J2/J1) to the 
multiply-degenerate state. How each bond type affects 
the gapped state established by the other is unknown, as 
both gaps may close, leaving some intermediate phase, or 
there may be a crossover between the gapped states. In 
the former case, our results provide a strong indication 
against the possibility that the gapless phase is ordered, 
but not that of a disordered, gapless phase in 2d arising 
due to the competition of the couplings. In the latter 
case, one may postulate a crossover where the singlet- 
triplet gap changes smoothly, remaining finite while the 
low-lying singlet states of the zig-zag chain limit are sim- 



ply split by increasing J 2 /Ji, or may also undergo a 
level-crossing with a higher- lying singlet. Examination 
of ground-state quantum numbers by ED will be partic- 
ularly useful in resolving this issue. 
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FIG. 5. Schematic phase diagram of the frustrated coupled 
ladder system. 

Fig. 5 summarizes the preceding analyses and discus- 
sion. The shaded regions represent those where open 
questions remain to be addressed. Returning to the ma- 
terials which have the frustrated coupled ladder struc- 
ture, the original ladder compound SrCu203 would ap- 
pear close to the isotropic point on the axis J\ = 0, and 
is well described as a liquid of resonating singlets residing 
primarily on the ladder rungs. That the frustrating inter- 
action is actually FM is of little consequence here. The 
zig-zag chain compound SrCu02 has very strong J2, and 
appears close to the upper vertex in Fig. 5; while J' 2 is ex- 
tremely small in this material, an interchain coupling ex- 
ponentially small in J2/J1 is sufficient to destroy the MG- 
type spin-gap state, so the real material may be located in 
the crossover regime of the previous paragraph. CaV205 
shows a large spin gap, a result in qualitative agreement 
with expectation for a system with J\ ~ J2 ~ J' 2 f-\ as we 
have seen that the frustrating interaction J\ must be sig- 
nificantly larger than the ladder terms in order to remove 
the system from the non-degenerate gapped phase. 

In summary, the frustrated coupled ladder system rep- 
resents well the wealth of interesting physics to be found 
in low-dimensional quantum magnets, and provides a 
valuable means of studying quantum phase transitions. 

We are indebted to M. Troyer for the quoted QMC 
result, and acknowledge also useful discussions with J. 
Bonvoisin, P. Millet and D. Poilblanc. 
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